Entanglement Skyrmions in multicomponent quantum Hall systems 
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We discuss charged topological spin textures in quantum Hall ferromagnets in which the electrons carry an 
pseudospin as well as the usual spin degree of freedom, as is the case in bilayer GaAs or monolayer graphene 
samples. We develop a theory which treats spin and pseudospin on a manifestly equal footing, which may also 
be of help in visualizing the relevant spin textures. We in particular consider the entanglement of spin and 
pseudospin in the presence of realistic anisotropies. An entanglement operator is introduced which generates 
families of degenerate Skyrmions with differing entanglement properties. We propose a local characterization 
of the latter, and touch on the role entangled Skyrmions play in the nuclear relaxation time of quantum Hall 
ferromagnets. 

PACS numbers: 73.43.Lp, 73.21.-b, 81.05.Uw 
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I. INTRODUCTION 

Quantum Hall systems with supplementary degrees of free- 
dom are fascinating objects of study, exhibiting phenomena 
such as the superfluid behaviour in quantum Hall bilayers at 
total filling factor v = iJ£iMiMd& In those systems, the su- 
perfluid ground state is isomorphous to a quantum Hall mono- 
layer ferromagnetic ground state at the same filling factor, but 
in terms of the 5=1/2 pseudospin degree of freedom en- 
coding which layer an electron finds itself in. A similar pseu- 
dospin ferromagnet is surmised to occur in a graphene mono- 
layer, where the pseudospin degree of freedom is analogously 
connected to the two degenerate mid band valleys i 9 ' 10 ' 11 

In quantum Hall ferromagnets at v = 1, the lowest charged 
excited states are topological defects, called Skyrmions, 12 
which in turn are fascinating in their own right. Indeed they 
are spin textures carrying a topological quantum number and 
which necessarily carry a quantized net charge. The proper- 
ties of Skyrmions in quantum Hall (QH) physics have been 
studied extensively since their theoretical prediction 13 and ex- 
perimental confirmation 14 in the mid 1990s. These charged 
topological spin textures are due to the strong electronic cor- 
relations once the low-energy excitations are restricted to a 
single Landau level (LL). They can exist as quasiparticles of 
a ferromagnetic ground state near a LL filling factor v = 1, 
where v — n e i/riB is the ratio between electronic density 
n e i and ns = eB/h = l/2irl B that of the flux quanta h/e 
threading an area of the two-dimensional (2D) electron gas 
at a particular magnetic field B (here, I b = yh/eB is the 
magnetic length). 

In bilayer QH systems, in contrast to the intrinsic electronic 
spin 1 /2 in monolayer QH systems, the effective interaction 
becomes pseudospin-dependent due to the different interac- 
tion strength between electrons in the same layer and that be- 
tween particles in different layers. It is indeed energeti- 
cally unfavorable to have all electrons in the same layer be- 
cause of a finite amount of charging energy per particle. This 
results in an easy-plane ferromagnetic ground state, i.e. in 
which the pseudospin magnetization is constrained in the xy- 
plane. Bilayer ferromagnets have, thus, an easy-plane [U(l)] 



symmetry. 

Topological excitations of easy-plane ferromagnets in bi- 
layer QH systems are conventionally described in terms of 
merons classified by their charge (±1/2) and their vorticity 
(il)^ 5 - As each meron has a logarithmically divergent energy, 
one needs two merons with opposite vorticity to make a topo- 
logically stable excitation of finite energy (bimeron), which 
may also be viewed as a pseudospin Skyrmion. 

Besides the physical spin in monolayer and the layer pseu- 
dospin in bilayer QH systems, valley indices can be encoded 
by pseudospin degrees of freedom. The indirect-gap semi- 
conductors Si and AlAs are well-studied instance,^i2£ as is 
the pair of Dirac points in graphene. 21 

One particularly interesting situation arises when the inter- 
nal Hilbert space comprises both the spin and the additional 
pseudospin. This combination can in principle lead to an in- 
ternal space with a high symmetry, such as SU(4) for a combi- 
nation of a spin and pseudospin pair of SU(2)'s. The proper- 
ties of SU(4) Skyrmions were considered by Ezawa 22 and that 
of general SU(iV) Skyrmions for > 2 by Arovas et alQ 
SU(4) Skyrmions have recently been revisited by Yang et al. 
in the context of graphene.— 

However, the Hamiltonian generically contains anisotropy 
terms violating the full symmetry. In the simplest situation 
of a Zeeman field, the symmetry o z <-* — o z is removed. 
Such symmetry-breaking terms can in principle also be ma- 
nipulated, for example the effective g-factor determining the 
strength of the Zeeman coupling can be tuned by applying ex- 
ternal pressure thanks to spin-orbit coupling. The energetics 
of the resulting topological excitations has been an active field 
of study in recent years . 25 ' 26 ' 27 ' 28 ' 29 ' 30 ' 31 

In this paper, we consider the properties of SU(2) x SU(2) 
Skyrmions in the presence of realistic anisotropy terms. We 
do this from a perspective of entanglement between the spin 
and pseudospin degrees of freedom. As our point of departure, 
we parametrise the general SU(4) Skyrmion using a Schmidt 
decomposition. 32 Our parametrisation is manifestly symmet- 
ric between the two SU(2) copies arising from spin and pseu- 
dospin. In addition, a third ('entanglement') spinor appears. 
These features turn out to be of help in visualizing the nature 
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of different types of SU(4) Skyrmions, for instance by pro- 
viding a simple picture for the relative entanglement between 
spin and pseudospin. 

With this in hand, we compute the Berry connection and 
then display in a transparent way the properties of entan- 
glement textures. (The topological stability of the SU(4) 
Skyrmion becomes, strictly speaking, apparent only in the as- 
sociated CP 3 description.) Crucially, we show that there exist 
entire U(l) families of degenerate Skyrmions which differ by 
their degree of entanglement. 

This in particular leads to the distinction of two different 
topological excitations, both of which have zero expectation 
value of transverse spin but in only one of which the spin de- 
gree of freedom is involved. At some special "maximally en- 
tangled" points, the 'order parameters' (a) and (f ) may even 
vanish altogether. 

This treatment offers some further insights into the mag- 
netic properties of multi-component systems, which have 
been studied experimentally in recent years, in the framework 
of nuclear magnetic resonance (NMR) experiments on bilayer 
QH systems^i 3 ^ In particular, we briefly discuss how entan- 
gled excitations affect the spin degrees of freedom measured 
in NMR experiments. 

The paper is organized as follows. In Sec. II, we intro- 
duce the basic model and the parameterization of the spin- 
pseudospin degrees of freedom. Spin-pseudospin entangle- 
ment is studied in Sec. Ill, where we propose a measure of 
such entanglement, the properties of which are described in 
Sec. IV. Sec. V reviews the relevant energy scales of bilayer 
and graphene QH systems, where entanglement Skyrmions 
may be physically relevant. The dynamical properties of en- 
tangled spin-pseudospin textures and their consequences for 
NMR measurements are discussed in Sec. VI. 



II. MODEL AND PARAMETERIZATION 

We study a two-dimensional electron system (2DES) in 
the QH regime, where the kinetic energy is quenched. We 
consider two internal degrees of freedom which provide the 
multicomponent nature of the 2DES. Firstly, the spin de- 
gree of freedom, represented by the Pauli matrices a a , where 
a = x,y, z. Secondly, the pseudospin, represented by a sec- 
ond set of Pauli matrices, r^, with /i = x,y,z. For con- 
creteness, it may help to have the bilayer index in mind when 
considering the pseudospin. 



A. Parameterization 

One may typically describe a texture at overall filling factor 
v = 1 by a Hartree-Fock state of the form 

|d>) = H[w 1 (X)c+ ht + w 2 (X)c+ tb + w 3 (X)c+ lt 
x 



+w A (X)c+ Ub ] |0) (1) 
where X labels orbital states in the plane restricted to the 
lowest LL. Internal states are labeled according to the projec- 
tion of the spin along the external magnetic field (f or J,) and 
the layer (i or b). We assume that the orbital states are each 
fairly well localized (on a typical scale given by the magnetic 
length) and mutually orthogonal. The explicit construction of 
such a basis for any Landau level is given for instance in Ref. 
|40| and has been used to describe the bilayer system by a lat- 
tice SU(4) model. 41 With such a choice, X = r can be viewed 
as the the location of the center of a given Wannier orbital on a 
von Neumann lattice. We will always assume that the texture 
varies slowly on the scale of the magnetic length, so that the 
Wi (1 < i < 4) may be regarded as smooth complex func- 
tions of a continuous position variable. We will also assume 
these local four-component spinors to be normalized, that is 
Yli—i \ w i( r )\ 2 = 1- To get a more direct physical interpreta- 
tion of such texture, it is rather natural to express the internal 
state at site r in a way that treats both spin and pseudospin 
degrees of freedom on an equal footing. This is achieved by 
the Schmidt decomposition 32 



|*(r)) = cos | |^s) !</</) 

+ sm|e^|x5)lx/> 



(2) 



where a and (3 are functions of r, and the local two- 
component spinors \if>s)> \xs), IV 1 /)* an d \xi) are constructed 
according to 



(3) 
(4) 

Here, 9 and <f> are the usual polar angles defining a vector 
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(sin# cos (j), sin 9 sin 0, cos ( 



(5) 



We will denote the pair for spin (pseudospin) by 9s,(f>s 
(6 1, 4>i), respectively. Note that the above parameterization 
for the two component spinors is not unique, because for a 
given classical unit vector n(9, <fi) we can multiply both \ip) 
and \x) by global phases that may depend on r. We have cho- 
sen the parameterization in order to minimize the occurrence 
of phase singularities - they appear only when one of the an- 
gles 9s, 9 j, or (3 is equal to tt. 

With these notations, the general four-component local 
spinor (Eq.|2]i has the form 
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|*(r)) 











w 3 




\w 4 / 





~T 2 e 



sin%sin^-sinfe l(/3 -* s " 

cos \ sin ^ sin ^ e i( ^~^ 
cos ^M- cos ^ sin 



(6) 



We thus have six parameters (63, 4>Si ®ii 4>ii a i P)> as ex ~ 
pected for four complex components minus a global phase and 
the overall normalization. All of these can vary as a function 
of the spatial coordinate r = (a;, y) in the xy-plane. 

The great advantage of factorizing the wavefunction in such 
a manner is that we can read off directly the reduced density 
matrices for the spin and pseudospin sectors: 

p s = ^(|V)^|) = cos 2 ||^)(^| + S in 2 || X s)(X5|, 

Pl = Tr s (|^)(V;|)=cos 2 ||^)(^|+sm 2 || X/ )( X/ |, 

for the spin and the pseudospin, respectively. This yields 

m a s = Tv{p s S a ) = cos a(i> s \S a \i>s) = cosan a {9 s ^s) 

(7) 

for the local spin and 

m£ = Tr(p/P^) = cos aiipilP^ipi) = cos a n»(6i, 0/) 

(8) 

for the local pseudospin density. Notice that for the case a ^ 
or 7r (i.e. cos 2 a < 1), the local (pseudo)spin densities are 
no longer normalized, but are of length |rri5/j| 2 = cos 2 a. 
Thus, in a semiclassical picture, the (pseudo)spin dynamics 
is no longer restricted to the surface of the Bloch sphere, but 
explores the entire volume enclosed by the sphere (Fig. Q]). 

B. Topological densities 

For an SU(2) symmetry, spin textures in the 2D plane may 
be classified according to which homotopy class they belong. 
First, one needs to map the xy plane, by stereographic pro- 
jection, onto the S 2 sphere. For an XYZ ferromagnet, the 
most general spin texture parameterization maps on the sur- 
face of the Bloch sphere S 2 . The mapings S 2 — ► S 2 , which 
can be continuously deformed into one another, may be col- 
lected into homotopy classes. 33 The latter form a homotopy 
group iT2(S 2 ) = 7L. The elements of TL, which are inte- 
gers Q, characterize different topological sectors. Spin waves, 
for instance, belong to the topological sector with topological 
charge Q = 0. Skyrmions carry a non-zero topological charge 
Q 0. In the QH ferromagnet, a spin texture generates a 
charge density connected to the former through the so-called 
topological density. The topological density arises because 
an electron travelling in a spin texture acquires a Berry phase 
analogous to an Aharonov-Bohm phase. The total charge, ex- 
pressed in units of the electron charge, of a spin texture in a 
QH state with Laughlin filling factor v = 1/ (2m + 1) is given 
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FIG. 1: Bloch spheres for entangled spin-pseudospin systems. Bloch sphere 
for the spin (a), pseudospin (b), and a third type of spin representing the en- 
tanglement (c). In the case of spin-pseudospin entanglement (| cosa| 7^ 1), 
the (pseudo)spin-magnetizations explore the interior of their spheres, respec- 
tively (black arrows). 

Next, we transpose these ideas to the larger SU(4) spin- 
pseudospin symmetry. With the help of the Schmidt decom- 
position, one may calculate the charge density induced by a 
general spin-pseudospin texture. When the internal states are 
described by a single SU(2) spin, the local density is sim- 
ply given by the topological density of the associated director 
field, that is 

Ptop = ^n(0, 0) ■ [0,n(0, 0) x 0-n(0, </>)], (9) 

where e 1 ^ is the antisymmetric tensor, e xy = —e vx = 1 and 
f n _ q g enera ii ze this to the present case, it is sufficient 
to note that the spatial variation of |^(r)) due to the texture 
induces a Berry phase, which is encoded by the Berry connec- 
tion, 

A(r) = -(¥|W>. (10) 

1 

The local charge density carried by the texture is then the 2D 
curl B(r) — [d x A y (r) — d y A x {r)]l (2vr) of the fictitious vec- 
tor potential A(r). After some simple algebra, we get 

A(r) = sin 2 + 
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Os Oj 
cos a ( sin 2 —V(j>s + sin 2 -^^4>i ) > (H) 



the curl of which is, thus, 

B(r) = cosa{ / 9 top [n(6's,(?!'s)] +p top [n(6'/,0/)]} + 

/j top [n(a,^)]+ (12) 

sin 2 — p top [n(a, <j>s)] + sin 2 —p top [ri{a, 0/)] . 



III. SPIN-PSEUDOSPIN ENTANGLEMENT 
A. A measure of entanglement 

The above parameterization makes it clear that there are 
Skyrmions in which spin and pseudospin can be perfectly 
entangled. This is most simply the case where for instance 
\^s) = (I), \xs) = (?), = (1/V2) (\) and | X/ ) = 
(l/y/2) ( i ) for all r, and the Skyrmion is a texture in a(r) 
and /3(r) only [which can take the same form as 9s, 4>s m a 
pure spin Skyrmion. This form is explicited in Eq. (l20t be- 
low]. In this case, it is easy to see, e.g., that a x = throughout 
as this operator is off-diagonal in this basis but has no matrix 
elements between \i/js)\ipi) an d |xs)|x/)- At the same time, 
o z has the same profile as in a standard spin Skyrmion, be- 
cause its expectation value is simply given by cos a(r). 

The vanishing of the expectation value of the transverse 
components of the spin is thus a consequence of their being 
entangled with the pseudospin degree of freedom. To quantify 
this entanglement, we propose to use the following measure: 



S: =i-5>r = i-£< 



-p,\2 



(13) 



Eq. ([T3l is easily expressed from the Schmidt decomposition 
Eq. © as 



a = sm a 



(14) 



As a result, for the abovementioned entanglement texture, 
which may be visualized on a third Bloch sphere [Fig. Q2 C )L 
local entanglement is minimal (S = 0) when sin a = 0, that 
is a — or a = it, corresponding to the center of the texture 
or to the points at infinity. Entanglement is maximal (H = 1) 
when a = ir/2, that is at the a distance from the origin given 
by the overal size A of the texture. Note that at a maximally 
entangled point, the moduli of both spin and pseudospin van- 
ish, as may be seen from Eqs. (0 and (8). 

Alternatively, for the spinor defined in Eq. ©, 



H = 4|wiW4 — W2W3Y 



(15) 



For a factorizable state, in which the wavefunction can assume 
the form of an outer product between spin and pseudospin 
parts, I^a/j) = u\v^, one immediately finds S = 0. Note that 
conversely, a vanishing 3 implies that we have a factorizable 
state. For a completely entangled state, e.g. w\ = W4 = 
1/V%, u>2 = wz = 0, we have S = 1. 



B. Degenerate texture families in the presence of spin 
anisotropy 

In the presence of a Zeeman field, the SU(2) symmetry of 
the spin degree of freedom gets reduced to a U(l) symmetry 
of rotations around the z-axis (along which we take the ex- 
ternal field to be aligned). The generator of this symmetry 
is a z . We furthermore consider, for illustration reasons, that 
the SU(2) pseudospin symmetry is likewise broken to U(l), 
and we choose the z-axis such that the rotations generated by 
t z , the z-componentof the pseudospin, commute with the un- 
derlying Hamiltonian. We discuss, in Sec. V, some physical 
examples, in the context of bilayer quantum Hall systems and 
graphene, where this situation arises. 

Given that the underlying model Hamiltonian commutes 
with t z and with a z , it also commutes with their product, 
a z r z . This means that a textured state |T) belongs to a de- 
generate family of textures |T(7)) = cxp(ijT z a z )\T). The 
most interesting property of this latter symmetry operation is 
its ability to transform a factorizable state into an entangled 
one. Indeed, it acts on the four component spinors accord- 
ing to 11^4 — > wi^exp(ij) and W2.3 — > 1^2,3 exp(— 17). 
In the process, \w1W4 — u^u^l 2 — > lu^u^l 2 + |w2i«3| 2 — 
[W1W4W2W3 exp(4i7) + c.c], where the bar indicates com- 
plex conjugation. Therefore, as 7 varies, S varies between 
Sm,» with 



4(|-u;i m\ T \w2W3\f 



(16) 



It is interesting to specify how one may, thus, create a max- 
imally entangled state. This corresponds to S„ iaa; = 1, or 
equivalently to |u>i| = |w4| and |u>2| = \ws\, because one 
may rewrite 

Smaz = 1- [(K| + K|) 2 + (|W 2 | + |W 3 |) 2 ] X 

[(Kl - Kl) 2 + (Kl - Kl) 2 ] ■ (17) 

Note that creation of a maximally entangled state starting from 
a factorized one is possible when the four spinor components 
have the same modulus. Another limiting case is obtained 
when Smin — ^max, that is the action of the <j z t z gener- 
ator does not change the degree of entanglement, although 
the quantum state does change as the angle 7 varies. This 
requires that at least one of the four spinor components van- 
ishes. If the initial state is factorizable, this translates into 
W1W4 = W2W3 = 0. 



IV. CONSEQUENCES OF ENTANGLEMENT 

A. Descendants of Skyrmions and bi-merons 

Let us now consider the effect of these entangling transfor- 
mations on the textures that are most relevant for the physics 
of the QH effect in a bilayer system at a total filling factor v 
close to 1. We choose, as a starting point, a 'bimeron' tex- 
ture, Ai, in which the spin is constant along the magnetic 
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field and the pseudospin has a topologically non-trivial tex- 
ture [see Eq.[28]below]. Due to the easy-plane anisotropy, the 
pseudospin is oriented in an arbitrary direction within the xy- 
plane far from the meron core, whereas in the core region it is 
oriented in the indirection. There are furthermore two pos- 
sible vorticities, such that one may distinguish between four 
meron types, each of which has half a charge. 15 A meron in 
itself costs a logarithmically large amount of energy due to 
the vorticity and must therefore be accompanied by a second 
meron of opposite vorticity, in order to obtain a pseudospin 
texture of finite energy. If both merons have the same topo- 
logical charge, one obtains a bimeron with an overall topo- 
logical charge of ±1, otherwise the bimeron is topologically 
connected to the (ferromagnetic) ground state. The action of 
t z on this state is simply to rotate the texture in the easy plane 
of the pseudospin, which is a global symmetry operation. The 
action of a z is trivial, as all spins point in the z-direction and 
hence a rotation about this axis has no effect. Therefore, the 
combined effect of t z <j z is the same global rotation as that 
effected by t z alone, and no entanglement is generated. This 
is in agreement with the general discussion given above, be- 
cause we have for instance W3 = W4 = if all the spins are 
up, which clearly satisfies W1W4 = W2W3 = 0. 

Next, we consider Sondhi et aJ.'s original spin-Skyrmion 
texture |<S), in which the pseudospin lies in a fixed direction 
in the xy-plane (the x-direction, say). The effect of a rotation 
generated by t z is simply to rotate the pseudospin uniformly 
about the z-axis in pseudospin space. The effect of a z rota- 
tions is a global rotation of the spin texture about the z-axis 
in spin space. Now, the crucial observation is that the effect 
of the combination t z o z is no longer just a trivial rotation. 
Rather, the operation generates a texture in which spin and 
pseudospin are entangled. To be more specific, we choose the 
pure spin texture as 



0(r) 



l«S> = 



sin ^rp-e* 



(r) 



Mi 



(18) 



where r = (x,y) denotes the location in the plane to which 
the local spinor \S) is associated. The second term in Eq. 
( TT8l precisely means that the pseudospin is homogeneously 
oriented in the x-direction, whereas the first one describes the 
(pure) spin texture. For a usual Skyrmion located at the origin, 
we have 



which shows that the complex number (x 4 iy)/\ is simply 
given by the stereographic projection from the south pole of 
the unit vector (er) onto the complex plane. 

Let us now apply the symmetry unitary operation 
exp[if (a 2 + 1)(t z - 1)] to \S). We, thus, obtain the state 
\£), which may be written as 



■ 0(r) i 

sm T e 



1 



(23) 



In this case, one recovers the pure entanglement texture men- 
tioned at the beginning of Sec. [Til] As we have already men- 
tioned there, (a x,y ) = in this state whereas the value of a z 
changes just as for an ordinary S spin texture. Likewise, we 
find that (t v > z ) = 0, and thus 



cos#(r) 



A 2 



A 2 



(24) 



B. General CP 3 Skyrmions 



Finally, it is interesting to extend this analysis to the most 
general CP 3 texture, subject to the condition that r x = — a z = 



1 for the 'point' 
\Q) defined by 



00. It is generated by the local spinor 



|a> = (|A 1 | 2 + |A 2 | 2 + 2|z| 2 + 2|6| 2 )- 




(25) 



where z = x 4 iy, and where we denote the normalization 
prefactor as JV in what follows. Here Ai, A2 and b are com- 
plex parameters. The quantity \b\ 2 4 (| A x | 2 + | A 2 | 2 )/2 is the 
overall length scale of the texture. The physical meaning of 
these parameters is seen more explicitly, by considering spe- 
cial cases, which allow one to recover the two textures (fT8l 
and (l23l already discussed, 



\S) = (2A 2 + 2|z| 2 )- 1 / 2 




(26) 



tan ■ 



(r) 



x + iy 
A 



(19) 



where A has the dimension of a length and corresponds to the 
spatial extension of the Skyrmion. In this fully separable state, 
the expectation value of the spin operator reads: 



2Aa 



A 2 


+ x 2 4 


V 2 




2Xy 




A 2 


4 x 2 4 


yl 


A 2 


2 

— x z — 


v 2 



A 2 



(20) 
(21) 
(22) 



\£) = (2A 2 4 2|z| 2 )- 1 / 2 



Ai 
-Ai 



and the abovementioned bimeron texture 



\M) = {2b 2 




(27) 



(28) 



Further insight may be obtained by considering the effect of 
a uniform rotation by an angle 7 of the spin around the z- 
axis. This changes (Ai,A2,&) into (e~ n \\,e~ 11 \2,b), after 
multiplying the spinor by a global phase that brings it back to 
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the above form. Similarly, a uniform rotation along the a;-axis 
in pseudospin space can be represented by the matrix 



ENERGY SCALES AND ANISOTROPIES IN BILAYERS 
QH SYSTEMS AND GRAPHENE 







e 2 cos ■ 




(29) 



acting on the column vector (Ai, A2, b) T . 

What are the entanglement properties of Ql We define 

2A = Ai + A2 and 26 = Ai — A2, so that we obtain 



[\5 2 z 2 \ + \X 2 b 2 \ + 2\X 6 zb\ cos cf)] 



(30) 



TV" 1 

where (j> — arg(z bXS). 

Thence, there are always two unentangled points, 5 = 0, 
the one by construction at\z\ = 00 as well as one at z m = 
—bX/5, at which the spinor is given by 



1 

lb 
Ai-A 2 



To maximize the entanglement, let us first consider \ z\ to be 
fixed. The above expression shows that we should pick <f> = 0, 
which means that z should belong to the line passing through 
z m and the origin, in the direction opposite to z m . It turns out 
that if z moves away from the origin along this half-line, 3 
first increases, reaches a maximum for \z\ = \zm\ and then 
decreases. \zm\ is given by: 



Z M 



JA6[ 



V¥F+W r +WJ+WWJW ,(3D 



which moves to 00 for 6 = Ai — A2 = 0. 

From this expression, one can read off the following special 
cases. The only completely unentangled textures, = 0, 
are those with (5 = and either b = or A = 0. The former 
corresponds to the Skyrmion S; the latter to the bimeron Ai. 

Complete entanglement at a point, E(zm) = 1, is achieved 
only in the case of |A| = \b\, for any 6. The case b = cor- 
responds to the entanglement Skyrmion £ introduced above. 
The special feature of b = is that maximal entanglement 
is obtained not just at two isolated points but on the circle 
I zm I = I A 1 1 . In that sense, £ is the most entangled texture. 

We may, furthermore, study the 'entanglement operator' 



r(7) = exp[i7T 2 (l + cr z )} , 
which transforms Q into Q' by replacing 

\[ = Ai exp[2i7] , X' 2 — X2 exp[— 217] 



(32) 



(33) 



Otherwise, the structure of the relevant equations remains un- 
changed, and A and S become 

A' = Acos(27) + i6sm(2-f) 
5' = 5cos(2-f) +iXsm(2-f) . 

These equations show that A' moves along an elliptical orbit 
in the complex plane. In the case where the modulus of b 
lies between those of the major and minor axes of the ellipse 
(A and S if they are real), there will thus be a value of 7 for 
which maximal entanglement occurs. 



A. Bilayer Hamiltonian 

In terms of the density operators defined in App. [A] the 
Hamiltonian for a bilayer quantum Hall system, taking into 
account both spin and pseudospin degrees of freedom, reads 

q q 

+A z S z {q = 0) + A T P x (q = 0) + A/P z (q = 0), 

(34) 



where the first two terms are due to the different intra- 
interlayer electron-electron interactions, with 



and 



-9 2 /2 



lie" 



eg 



in terms of the layer separation d and the dielectric constant 
e. The last three terms in Eq. d34i > are Zeeman-type terms; 
the first one is indeed the Zeeman term that acts on the z- 
component of the spin density, whereas the second one is due 
to interlayer tunneling, and the last one is a layer imbalance 
term which may be tuned by applying a gate voltage. Whereas 
the Zeeman-type terms break the SU(4) symmetry explicitly, 
the first interaction term is SU(4)-symmetric, and the second 
one gives rise to an easy-plane anisotropy for the pseudospin 
magnetization, and thus breaks the SU(2) pseudospin symme- 
try down to U(l). 

In the case of GaAs heterostructures, we have the energy 
scales (for a dielectric constant of e ~ 13 and an effective 
spin coupling g = —0.44) 



energy 


value for arbitrary B 


value for B — 6T (u = 1) 


e 2 /el B 




122K 


A z 


0.33B[T]K 


2K 


A T 




0... 100K 



where the last term is sample dependent. It can, for instance, 
be varied by considering samples with different distance be- 
tween the layers. Indeed, the tunneling term At can vary 
over a rather large range and - most importantly - it may 
become the smallest of these energy scales. One notices 
that the largest energy scale is given by the interaction en- 
ergy. Strictly speaking, the characteristic intralayer interac- 
tion is given by e 2 /els, whereas the interlayer correlations 
are governed by e 2 /e^d 2 + l 2 B . However, we typically have 

d ~ Ib = 26/ y/ B[T] nm. The Zeeman term is roughly two 
orders of magnitude smaller than the characteristic interaction 
energy. This means that although both the Zeeman effect and 
interlayer tunneling tend to fully polarize the system, with the 
spin aligned along the z-axis and the pseudospin along the x- 
axis, this polarization only occurs subsequent to the onset of 
ferromagnetic order. 
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Notice that the Hamiltonian ( [34-b also applies to the case 
of graphene, a system which is, however, less anisotropic than 
bilayerGaAs heterostructures. The interactions are dominated 
by the symmetric part, V+(q), of the interaction, whereas 
symmetry-breaking interactions are suppressed by a factor 
a/Is, in terms of the carbon-carbon distance a = 0.14 
nm / ' 11 ' 36 which replaces d as the second characteristic length 
scale. As for GaAs heterostructures, the Zeeman effect in 



graphene plays a minor role as compared to the leading inter- 
action energy scale. Concerning possible Zeeman-type terms 
acting on the pseudospin, which describes the two inequiva- 
lent Dirac points, their existence remains an open and intense 
field of research i 37 ' 38 i 39 However, these terms are expected to 
be also small with respect to the interaction energy. The en- 
ergy scales for graphene are summarized in the following tab- 
ular (for typical values e ~ 2.5 and g ~ 2) 



J 



energy 


value for arbitrary B 


value for B — 6T 


value for B = 25 T 


e 2 /d B 




620K 


1250K 


A z 


1.2S[T]K 


7K 


30K 


A sb < {e z /d B )(a/l B ) 


< 1B[T]K 


< 6K 


< 25K 



Energy scales [K] (for bilayers at 6T) 



Energy scales [K] (for graphene at 25T) 



47tp s 



o- 



-30 



~2 
< 1 



T z , G a , T z <7 a 



T z , <7 Z , T Z <7 Z 



FIG. 2: Energy scales for GaAs bilayer QH system for a magnetic field of 
B = 6 T. The spin stiffness is given by pg = (16v / 27r) _1 e 2 /(/Be). On the 
RHS, we indicate the symmetry generators for the Hamiltonian obtained by 
neglecting terms with energy scales below fcg T. 



We have displayed the values both for 6T and for 25T because 
the filling factor v = 1 may be obtained for different values 
of the magnetic field, due to a control of the carrier density 
by application of a gate voltage. The latter choice (25T) is 
motivated by the fact that above this value, the v — 1 QH 
plateau is well developed.— 



B. Physical relevance of entanglement in bilayer systems and 
graphene 

The above considerations and the possibility of entangled 
texture states sensitively depend on the smallness of the tun- 
neling gap in bilayer systems. Experimentally, a small tun- 
neling gap is a reasonable assumption as discussed above, and 
we consider the case At = in the following. 

In the temperature regime T < Az, one may then find 



47ip s 



• sb 







380 



-30 
<25 



full SU(4) symmetry 
%^ o a x^o a 



h\ ? 



FIG. 3: Energy scales for graphene in B = 25 T. On the RHS, we indicate 
the symmetry generators for the Hamiltonian obtained by neglecting terms 
with energy scales below kgT. 



entangled states with an internal symmetry spanned by the 
Cartan algebra {t z , a z , r z a z } of the SU(4) group. As shown 
above, the degeneracy supports entangled texture states in this 
regime. In the absence of Az, further entangled states arise, 
described by an SU(2) x U(l) symmetry group, which is 
larger and non-Abelian - the operators {t z , a a , T z a a }, do not 
all commute due to the algebra satisfied by the Pauli matrices 
a a (see Appendix A). This means that one may chose any a a , 
instead of just a z , to create entangled states, with the help of 
the entanglement operator 

r( 7 a )=exp[z 7 V z (l + a a )], 

which generalizes that in Eq. d32l . A full analysis of 
Skyrmions in this regime is deferred to future work. 

The graphene case is quite similar. The energy scales are 
shown in Fig. [3] for a characteristic field of B = 25 T, where 
the v — 1 plateau in the Hall resistance is well established. 49 
The interaction potential has been shown to be SU(4) symmet- 
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ric over a large temperature range, and a broken SU(2) pseu- 
dospin symmetry due to lattice effects is expected only below 
an energy scale A sb ~ (e 2 / 'd B ){a/l B ) ~ B[T] K . 10 i"- 36 i 38 
The nature of this symmetry breaking differs in the relativis- 
tic LLs | n | > 1 from that in the central one n = 0. For 
\n\ > 1, one expects an easy-plane anisotropy, 10 which sup- 
ports entangled states with an internal symmetry generated 
by {t z , a z , r z a z }, similarly to the bilayer case in the regime 
At < T < Az- The n = case is more involved, and the 
easy-plane anisotropy of the pseudospin magnetizationiSi 3 - 6 . is 
in competition with an easy-axis anisotropy, 11 ' 38 which may 
lead to a residual Z2 pseudospin symmetry. However, in all 
cases numerical prefactors strongly suppress the energy scale 
for lattice effects A s b ~ (e 2 /e/s)(a/?s) below the typical 
Zeeman energy scale Az ~ 1.2S[T] K. 

In an intermediate regime of A s & < T < Az, the in- 
ternal symmetry is given by {t^, a z , T^a z }, which supports 
entangled texture states. For temperatures well above Az, 
one recovers in turn the full SU(4) symmetries, and SU(4) 
Skyrmions 22 ' 23 ' 24 are expected to be the relevant charged ex- 
citations. 



VI. DYNAMICAL PROPERTIES OF ENTANGLED 
TEXTURES 

One strong experimental motivation for the present work 
has been the experimental observation of nuclear spin relax- 
ation 34,35 in the vicinity of the transition between the quan- 
tum Hall pseudospin ferromagnetic phase and the compress- 
ible phase in a bilayer system at total filling factor v = 1, 
where the control parameter is the ratio d/ls ■ These two stud- 
ies agree on the fact that nuclear relaxation is much faster in 
the compressible phase than in the quantum Hall phase. What 
is particularly striking is that a slight deviation of the total 
filling factor away from v = 1 in the quantum Hall phase in- 
duces a significant increase in the relaxation rate \jT\. Such 
a strong variation with the filling factor is reminiscent of what 
has been observed in a monolayer at v — l, 14 where the ad- 
ditional relaxation is attributed to the gapless spin modes in- 
duced by the presence of Skyrmions. 26 In the context of a bi- 
layer, this was surprising, because the samples are widely be- 
lieved to be fully spin-polarized in their quantum Hall phase 
at v = 1, and therefore the charged textures would be pseu- 
dospin bimerons, which do not couple to nuclear spins. The 
dominant interpretation of these data is based on the idea of 
a coexistence between both quantum Hall and compressible 
phases in the vicinity of the transitions 4 ^ This idea has been 
first proposed^ 2 - to explain the strong longitudinal Coulomb 
drag peak observed in the vicinity of the transition 18 and is 
supported, on the theoretical side, by numerical diagonaliza- 
tion studies 43 which point towards a first order transition. Fur- 
thermore, some recent experiments have confirmed the fact 
that there exists a large jump in the spin polarization between 
the two phases^ 4 . 

Here we would like to revive the case for another scenario. 
It is not clear to us why the compressible phases found either 
at v = 1 and large d/ls or away from v = 1 for smaller 



d/ls should have the same spin excitations. This assumption 
is supported by experimental results 34,35 which show a com- 
parable 1/Ti in both cases (see for instance figures 3(b) and 
3(c) in Ref. 34). But theoretically, we may question why com- 
posite fermion physics 45 which is expected to emerge at large 
d/ls with each layer at half filling, should remain relevant at 
small d/ls and away from v — 1/2 per layer. Note however 
that numerica l 46,47 and analytical 48 calculations indicate that, 
even for a homogeneous system, a small but non-zero layer 
separation admixes composite-fermion components to the in- 
compressible superfluid state. An alternative way to account 
for the strong relaxation away from v = 1 in the gapped phase 
is to assume the presence of entangled textures. On purely en- 
ergetic grounds, the pure spin texture is not favorable, because 
of the large Zeeman energy. We may then expect that an en- 
tangled texture would have a significantly lower energy than 
a spin texture, and that it can still contribute to nuclear spin 
relaxation. It has in fact been shown that a lattice of entan- 
gled skyrmions has even a lower energy than a lattice of spin- 
polarized bimerons if d/ls is large enough and the Zeeman 
coupling is not too large. 30 Leaving aside these delicate ener- 
getic questions, which require to optimize the CP 3 profile in 
the presence of the various anisotropies, 29,30 we now briefly 
address the question of the collective excitations around an 
entangled texture, and their impact on electronic spin-spin 
correlation fucntions. Note that detailed numerical calcula- 
tions of the collective modes for CP 3 skyrmion crystals are 
available^ Here, we would like simply to present a simple 
framework to analyze the small fluctuations around an entan- 
gled texture, focussing on the counting of zero modes. A strik- 
ing result is that the NMR relaxation rate induced by spin fluc- 
tuations around a given texture is changed when this texture 
is modified by the application of the entanglement operator 
exp(i"fa z T z ), even when this operator is an exact symmetry 
of the Hamiltonian, i.e. when [see Eq. (l34l >l vanishes. 

For simplicity, we consider the SU(4) symmetric case. In 
this situation, the general texture that carries a unit topological 
charge has the formSi 

w(z) — zu + v, (35) 

where w(z) parameterizes in fact the ray in CP 3 that con- 
tains this vector. Because both u = (ui, 7*2, 113, 114) and 
v = (t>i, t>2, V3, V4) have four complex components, we have a 
seven complex parameter family of degenerate textures, after 
factoring out the simultaneous multiplication of u and uby the 
same arbitrary complex number. We therefore expect seven 
zero modes for the dynamics of small fluctuations around such 
texture. Note that in the absence of the texture, w(z) is con- 
stant, so there are three zero modes that may be attributed to 
the ferromagnetic vacuum sector and four to the texture: one 
of them corresponds to a translation (z — > z + a) and the three 
remaining ones to the three independent directions in the in- 
ternal space of CP 3 . 

It is in fact illuminating to consider the case of q textures. 
Then, the minimal energy configurations are described as 50 

9 

3=0 
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where the iV-component q + 1 complex vectors Uj are arbi- 
trary. We have therefore a highly deformable system where 
(q + 1)N — 1 degrees of freedom exhibit no restoring force. 
For clarity, we are generalizing here to the CP W_1 model with 
global SU(A0 symmetry, N being equal to four in the case 
of a bilayer. Removing the N — 1 degrees of freedom as- 
sociated with the ferromagnetic vacuum sector leaves qN of 
them attached to the q textures, so finally, there are N static 
degrees of freedom per texture. Note that the SU(iV) symmet- 
ric CP^ -1 model describes a very special system, where the 
textures do not interact. In reality, various anisotropics and 
also the long-range Coulomb interaction between the elec- 
tric charges bound to the textures lift this massive degeneracy 
within this ideal g-skyrmion manifold, stabilizing most likely 
a crystal of skyrmions. Without anisotropies, nor Coulomb 
repulsion, the qN static degrees of freedom describing the 
small deformations around such crystal can be organized into 
TV non-dispersive Bloch bands. In the presence of realis- 
tic anisotropies and the Coulomb repulsion, these bands ac- 
quire a dispersive character. Clearly, one of these N branches 
of the collective spectrum will be a magnetic-field phonon 
mode, with 51 uj ~ |g| 3 / 2 , and N — 1 other branches will be 
Goldstone-type modes. We therefore recover the counting that 
emerged from a more microscopic analysis^ 

Let us now return to the quantum dynamics of the electron 
system around a single texture. As in Refs. 1, we 

consider the framework of the time-dependent Hartree-Fock 
(TDHF) approximation, which can be formulated in several 
different ways. A rather appealing one consists of viewing 
the set of Slater determinants as a classical phase spaced The 
Hartree-Fock Hamiltonian is then used as a classical Hamilto- 
nian that generates the same dynamics as the TDHF equations 
of motion. In our problem of a filled LL, we also assume 
(see the beginning of Sec. II) that the occupied single-particle 
states correspond to localized four-component spinors u)j (r) 
that are not normalized here. Indeed, dropping the normaliza- 
tion constraint simplifies the equation of motion at the expense 
of introducing a form of a generalized local gauge symmetry 
Wj(r) — > h(r)wj(r) which has to be factored out in the over- 
all mode counting. The Hartree-Fock Hamiltonian may be 
taken [assuming SU(4) symmetry] as 



H({w}) 



1 



6fr 



(Vw, Vto) (Vm,w)(w,Vw) 



(id, w) 



(w, w) 



(37) 

where the inner product is defined by (w,w') = 

Et=i ^j{ r ) w j( r )- Thi s i s me onr y l° ca l Hamiltonian that is 
quadratic in gradients, invariant under the global SU(4) trans- 
formations in spinor space and also under local gauge trans- 
formations Wj(r) — > h(r)wj(r)£& As usual, it is now conve- 
nient to introduce the complex position coordinate z = x + iy, 
such that 



H({w}) = = I d 2 r 



(dw,dw) + (dw,dw) 



(w, w) 

(dw, w)(w, dw) + (dw, w)(w, dw) 
(w, w) 2 

The notations mean in particular that (dw,w'] 



J2j(dz w j) w j an( l (dw,w') = Y^j(BzWj)w'j. The classical 
Hamilton equations of motions are 



d , . SH 
ot °Wj(r) 



■ d - ( 



5H 
8wj(r) 



(38) 



Because of the local gauge invariance, the norm (u>(r), w(r)) 
is conserved for any r. The price to pay for this is that the 
equations of motion are rather complicated. Using the gauge 
symmetry, we may modify Eq. ( 1381 ) by a term proportional 
to Wj (r) that simplifies it considerably, giving up the conser- 
vation of (w(r),w(r)), while preserving the ray of w(r) in 
CP 3 and therefore the physical content of the solutions. This 
simplified equation of motion reads 

9 *2 Aw,d s w) (w,d z w) 
i—w = —o r ^w+— —o z w+— —OzW+\(r 7 t)w , 



dt 



(w, w) 



(w, w) 



(39) 



where A(r, t) is an arbitrary complex function of r and t, and 
d 2 z = d z dg- Note that the gauge transformation w'(r, t) = 
h(r, i)w(r, t) sends a solution w into another solution w' pro- 
vided A is changed into A' given by 



A' = x + ir 1 

_ (w, d z w) 
(w, w) 



dh 2 {w,d- z w) 

% lu + d *,s h ~ ~( r dzh 

ot (W, w) 



d- z h 



- 2- 



dghd.h 



h 2 



It is immediate to check that the general static texture w(r) = 
zu + v is indeed a solution with A = 0. 

Linearizing the equations of motion around these static so- 
lutions gives 



(w, d z w) 
[w, w) 



o z e+ — — d z e + \{r,t)w , 

(uj, w) 

(40) 

where e(r, t) denotes an infinitesimal variation and we have 
used the fact that d z w = for the static texture. Note that A 
is an infinitesimal complex function of the same order as e. A 
first order infinitesimal gauge transformation with h(r,t) = 
1 + £(r, t) acts as 



e' = e + (w , 
8( 



(w, d z w) 
(w, w) 



The eigenmodes have the form e(r, t) = e(r)e~ luJt . Because 
the classical Hamiltonian is bounded below, u> can take only 
non-negative values. The zero modes are obtained from the 
condition d z e = 0, and the constraint of a fixed topological 
charge imposes that e should be a polynomial of degree 1 in 
z. This yields an eight-dimensional space of zero modes in 
which the pure gauge transformations e = (ow with £0 con- 
stant in space and time, have to be removed. We therefore 
recover the seven dimensional complex space of zero modes 
discussed at the beginning of this section. 
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Note that we have a conserved Hermitian form fl(r), £) for 
any pair rj, £ of solutions of the linearized equations of mo- 
tion (l40l . This form ft is defined by 



7 j 



07,0 



(w,w) (w,w) 2 



(41) 

and is invariant under any local gauge transformation that af- 
fects simultaneously w, r\, and £. In particular, for eigen- 
modes 77(1-, i) = ri(v)e~ %UJ,lt and £(r, t) = ^(r)e~ 1 " 5 *, we 
have f2(?7,£) = as soon as w,, ^ u>£. This shows that a 
complete eigenmode basis e Q '"(r) (a being an internal polar- 
ization label) can be chosen to be orthogonal for the Hermitian 
form f2. 

We will not try to find the full continuous spectrum of 
the scattering problem defined by Eq. d40b . Instead, we 
shall now consider the effect of the entanglement operator 
exp(i"fa z T z ) on the local spin correlation function C(r, u>) — 



fdt(S+(r,t)S-(r,t = Q))e iuJt that controls the local nu- 
clear spin relaxation rate at the position r. Here, S ± (r, t) — 
S x (r,t) ± iS v (r,t), in terms of the local spin densities 
S x l y (r,t). In the presence of a texture, translation invari- 
ance is lost, so this correlation function will explicitly de- 
pend on r. The TDHF approximation does not give us 
directly wave-functions for collective modes, so we use it 
to express the corresponding response function R(r, u>) = 
| / dt{[S+(t), S-{0)])e iuJt in terms of a complete basis of 
eigenmodes e a ' u . R(r,oj) is related to the correlation func- 
tion by the usual fluctuation-dissipation formula, 



2h 



1 - e-P hi 



(42) 



sponse function is then expressed as 



Introducing the 4x4 matrices M ± = <J^, in <8> li so > the re- 



(w(r),M+e a ^(r))(e a ^(r),M-w(r)) 
(w(r), w(r)) 2 



e(io) 



(w(r),M- 



J (r))(^-(r),Af+ W (r)) i 



(w(r), w(r)) 2 



_J43) 



Notice that we have treated, within the calculation of the 
low-energy modes, the Skyrmion as a classical object, which 
is a valid approximation for large Skyrmions, A ^> Ib- 
For smaller Skyrmions, the translational modes may still be 
treated on the classical level, whereas the rotational modes 
need to be quantized. 16 

Now, let us start from the pure spin texture discussed 
in Sec. IV and let us apply the entanglement operator 
ex-p(i-fcr z T z ) to it. This unitary transformation in the full 
quantum Hilbert space induces a canonical transformation on 
the classical phase-space of Slater determinants. It is, there- 
fore, equivalent to compute the correlation function (S + S~) 
around the transformed texture and the correlation func- 
tion (S + ( r y)S~( r y)) around the original one, where S(7) = 
exp(— ija z T z )S exp(i-fcr z T z ). Explicitly, we obtain 

S±( 7 ) = cos(27)S* ± <g) 1 qp i sin(2 7 ),S' ± <g) t z . (44) 

For the fluctuations around the pure spin texture, a little in- 
spection shows that the t x component of pseudospin is con- 
served. Therefore, crossed terms involving only one r 2 oper- 
ator do not contribute to the sum over eigenmodes d43b . We 
finally get 

{S+S-) 1 = cos 2 (2 7 )(S+S-) 7=0 

+ sm 2 (2 7 )((5+T 2 )(5-T 2 )) 7=0 • (45) 

Note that the two correlation functions that appear on the 
RHS of the above equation will in general differ, because we 
are considering local spin and pseudospin operators. So the 
two t z operators taken at two different times are not expected 



to commute, in spite of the fact that the z-component of the 
total pseudospin is conserved. This expression is quite inter- 
esting because it exhibits a competition between two effects. 
When 7 increases slightly away from zero, the first term de- 
creases, in agreement with physical intuition - the local spin 
is reduced by entanglement, and so is its contribution to NMR 
relaxation. Such scenario would be consistent with the experi- 
mental observation that 1 jT\ increases faster for a monolayer 
than for a bilayer as the filling factor is moved away from 
the v — 1 value (see for instance Fig. 4 in l35l) . But this 
is not the only effect induced by entanglement - the second 
term also builds up, showing that new relaxation channels oc- 
cur, that correspond to flipping simultaneously both spin and 
pseudospin degrees of freedom. 



A physical consequence of this variation of l/T\ with the 
entanglement parameter 7 could be the existence of large fluc- 
tuations in the nuclear relaxation rate from one sample to an- 
other, or for the same sample after thermal cycling between 
the low temperature quantum Hall state and the normal elec- 
tron fluid. Indeed, we may speculate that, in the vicinity of 
the bilayer phase transition, the energy landscape in the space 
of possible textures is characterized by a rather flat minimum 
around the bimeron texture, and that it may then be diffi- 
cult for the system to relax to its absolute ground-state. If 
this picture were to be confirmed, the above discussion shows 
that 1/Ti would depend on the actual degree of entanglement 
reached by the system. 
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VII. SUMMARY AND OUTLOOK 
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We have investigated the possible entanglement between 
two internal SU(2) degrees of freedom of electrons restricted 
to a single LL. Whereas the physical spin constitutes one 
SU(2) copy, the second one describes a pseudospin degree of 
freedom such as the layer index in bilayer QH systems or the 
valley degeneracy in graphene in a strong magnetic field. The 
Schmidt decomposition allows one to treat both on an equal 
footing, while embedding them in a larger SU(4) space. This 
has enabled us to study the effects of realistic anisotropies in 
a transparent manner. In particular, by defining an entangle- 
ment operator which generates entire families of degenerate 
Skyrmions, we are able to construct wavefunctions for entan- 
glement Skyrmions which are easy to visualize. 

It has also become apparent that the NMR response of such 
multicomponent systems is affected by the degree of entan- 
glement. Indeed, it remains to be seen whether this may pro- 
vide an alternative to the puddle model of coexistence of a 
spin-polarized incompressible and a partially polarized com- 
pressible metallic phase, 42 which is mainly invoked in the ex- 
planation of the increased NMR relaxation rate when driving 
the total filling factor away from vt — l.-^ 4 ^ Unfortunately, 
in graphene, where anisotropies in the valley sector are quite 
weak, NMR measurements turn out to be difficult because the 
majority (~ 99%) of the carbon atoms are 12 C isotopes, which 
are NMR inactive due to a zero nuclear spin. 

This work leaves many open questions. The first is to iden- 
tify precisely the regimes in which entangled textures can be 
stabilized as the absolute energy minima for charged excita- 
tions. This has been demonstrated to occur in the presence 
of charge imbalance between the two layers^ or even for 
the balanced system, for a layer separation lower than, but 
close to the critical one, provided the Zeeman energy is not 
too large. 30 An important issue here would be to understand 
the effect of random impurities that may disorder such lat- 
tices of entangled skyrmions and thus modify the energy bal- 
ance between various configurations. But even when the spin 
polarized bimeron is the most favorable texture, it would be 
very interesting to understand better the energy landscape in 
its vicinity, because available numerical results suggest it may 
be rather flat, 30 which raises the possibility of a slow dynam- 
ics and a lack of equilibration. Finally, there exist certainly 
many other physical signatures of spin-pseudospin entangle- 
ment besides the local reduction in the magnitude of the av- 
erage spin and pseudospin vectors. These issues are left for 
future studies. 



APPENDIX A: LOWEST-LL ALGEBRA FOR 
MULTICOMPONENT SYSTEMS 

The local spin and pseudospin densities may be expressed 
in terms of Pauli matrices as 



and 



S a (r) = I/(r) ® (<r° ® 1) 



p/*(r) = i/(r) ® (1 ® t"). 



(Al) 



(A2) 



Similarly, the total (charge) density may be written as 

p(r) = /(r) ® (1 ® 1). (A3) 

Here, /(r) is the U(l) one-particle density projected to the 
lowest LL, the Fourier components of which satisfy the 
magnetic-translation algebra 1 *^ 



and 



[/(q),/(q')] = 2isin 



{/(q),/(q')}=2cos 



qAq' 



qAq' 



/(q + q') (A4) 



/(q + q'). (A5) 



In Eqs. (IAlt - dA3l l, the bar indicates that the density oper- 
ators are restricted to the lowest LL. Because of Eqs. ( IA4b 
and ( IA5b , the Fourier-tranformed density operators satisfy the 
commutation relations 

[p(q), p(q')] = 2* sin p(q + q'), (A6) 



[S a (q), p(q')] = 2i sin ( 5 a (q + q'), (A7) 



[P"(q),p(q')]=2isin 



qAq' 



P"(q + q') (A8) 



[5 Q (q),5 b (q')] = ^ fc sin' q f) 



p(q + q') (A9) 
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and 



[P^q),P-(qO] = ^sm(^^p(q + q') (A10) 

cos | SAS^ P CT (q + q'). 



These commutation relations are directly obtained from the 
relation 

[/( q )®A/(q')®£] = j([/(q)J(q')]®M,B} 

+{/(q),/(q')}®IAB]), 

(All) 

for any pair of operators A, B. Furthermore, one obtains the 
mixed spin-pseudospin commutator, 

[S a (q), P>(q')] = l - sin (^^) /(q + q') ® K ® r") . 

(A12) 

Naively, one might have expected that the spin and pseudospin 
densities are decoupled because of [a a ®l, 1(g) r' 1 ] = 0. How- 
ever, the last commutation relation indicates that the spin and 



pseudospin dynamics are indeed entangled when taking into 
account local variations of the texture. This is due to the pro- 
jection into the lowest LL, which yields the noncommutativity 
of the orbital degrees of freedom [Eq. (IA4H . It may also be 
seen directly from Eq. (IA lib : if we had [/(q),/(q')] = 0, 
the anti-commutator term ~ {A, B} would vanish, and one 
would obtain [S' a (q), P M (q')] = 0. The same argument also 
yields the spin-charge and the pseudospin-charge entangle- 
ment, revealed by Eqs. iAH and ( IA8b . 



Eq. (IA12b indicates that we generate, via the commuta- 
tors, the generators a a ® r p which, together with a a (g) 1 and 
1 ® r M , give rise to the larger internal symmetry group SU(4), 
in which the SU(2)xSU(2) group thus needs to be embed- 
ded. The SU(JV) extension of the magnetic translation group, 
Tm x SU(-/V), with N > 4, was indeed the starting point 
of several other theoretical works on bilayer quantum Hall 



systems 



22.23.28.29 



We have thus shown that a SU(4) descrip- 



tion of the texture states is necessary even if the ground state 
has the reduced internal symmetry SU(2)xSU(2). 
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